The interval description of dynamics of celestial bodies in the planetary problem 
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The interval approach to computation of dynamics of celestial bodies in the planetary problem 
has been considered. It is based on the refusal from idealization of infinitely high resolving capacity 
of measuring tools, and forms an absolutely exact algorithm free of round-off error accumulation 
effect. The possibilities of the proposed approach are shown by the examples of Kepler's Problem 
and the problem of stability of the Solar system major planets for time interval of 6 billion years. 
The comparison of the interval and classical predictions of Kepler's particle location in Kepler's orbit 
provides support for the effect predicted by the theory, namely - conservation of the interval within 
which the values of difference of interval and classical coordinates lay with time. The computational 
results of the Solar system major planet orbital dynamics agree with the results obtained with the 
classical approach. 
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I. INTRODUCTION 

The interval-discrete concept of dynamics of point sys- 
tems has been formulated in the paper The study 
is based on the refusal of idealization of an infinitesi- 
mal error of observations and calculations that is on the 
property of the limited resolving capacity of measuring 
tools and measurement processing facilities, which is not 
considered obviously in classical (trajectory) concept of 
motion. [13fl In the paper [2| the interval description of 
macro bodies motion is applied to the planetary prob- 
lem. The numerical computations of the Solar system 
planetary orbit evolution for time intervals of about 500- 
million years are performed in the study with the help of 
interval equations. The computations have verified the 
results obtained by classical methods, and have enabled 
to make a conclusion that the proposed interval approach 
is applicable for solving the problem of many bodies. 

This paper presents the continuation of the study hav- 
ing been begun in The study objective is to com- 
pletely realize and show the advantages of the interval ap- 
proach as applied to the planetary problem and, thereby, 
to lay the groundwork for experimental check of the equa- 
tions of the interval theory. 

In this connection we remind P, that formally the 
interval equations are a system of integer mappings of 
recurrent-type obtained by a special procedure of quan- 
tization of the time and spatial continuums of the system 
and the intervals of its dynamic variables. The feature of 
these mappings lies in the fact that they form an abso- 
lutely exact computational algorithm free of the effect of 
the round-off error accumulation. In this case it is pos- 
sible to carry out computations with an unlimited num- 
ber of iterations. The other feature of these mappings is 
a reasonable simplicity significantly reducing amount of 
computations in comparison with the classical approach. 
In the paper the mentioned features are illustrated by 
the examples of Kepler's Problem and the problem of 
stability of the Solar system for cosmogonic times. 

Kepler's Problem is of interest as a basic problem. In 



particular, it allows demonstrating the interval tube ef- 
fect 0| . Its manifestation is that the values of deviations 
of variable interval centers from their classical values lies 
in a strictly fixed interval the width of which does not ex- 
ceed the initially preset width of interval variables. Clas- 
sically it means that precision of the interval prediction 
for Kepler's particles dynamics remains a constant value 
irrespective of "integration" duration. It is natural that 
the real reliability of such prediction will be limited by 
the time of system's fall outside the limits of its «hori- 
zon of predictability)) [3|. In this connection the interval 
dynamics in contrast to the classical dynamics in which 
there is no concept of "horizon of predictability)) at all, 
allows to exactly estimate the "lifetime" of theoretical 
prediction, which is always finite in practice. 

Kepler's Problem allows demonstrating one more im- 
portant property of the interval dynamics. That is abil- 
ity of the interval tube for closing if the corresponding 
classical trajectory represents a closed phase curve. The 
uniqueness of this property is manifested by its inhesion 
to particularly numerical result, i.e. it means an abso- 
lutely exact recurrence of numerical parameters of the 
tube through one or several turns of a particle. 

And, at last, there is one more interesting effect, which 
is demonstrated by the example of Kepler's Problem. An 
interval particle alongside with an orbital moment has its 
own kinetic moment (spin). This effect is one of conse- 
quences of the interval nature of physical system space. 
It is conditioned by the fact that an interval particle in 
contrast to a classical one has the status of not a mathe- 
matical point, but a physical one. In other words the par- 
ticle has finite size (determined by the width of intervals 
of the problem spatial variables) and, as a consequence, 
possesses its own moment of an impulse. 

Freedom of the effect of round-off error accumulation 
allows successful solving of particularly academic prob- 
lems alongside with practical ones. One of the problems 
is a study of a dynamic stability of the Solar system ma- 
jor planets. Difficulty of its solution is caused, first of 
all, by long-term character of the required prediction of 
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the planets motion assuming carrying out the computa- 
tions with time intervals of the Solar system age order. 
In the present work such computations are made for a 
time interval of six billion years. The obtained results 
agree with the results obtained earlier within the limits 
of the classical approach as for external planets 0-01 , so 
as for internal planets [1, [H, 0] ■ 



integer centers r n ,cp n and Pn^ot the intervals (2) can be 
written down in the form of 



where 



r n +i = r n + [wp n 

Pn+l = [Pn+l], 
Pn+l = ['Pn+l], 



(4) 



II. KEPLER'S PROBLEM 

Let's consider an interval statement of Kepler's Prob- 
lem in the space of polar coordinates presenting Hamil- 
tonian of a classical analogue of the system under con- 
sideration in the following form: 



,■2 



H = l(Pr + ) - --{P V = c = const) 



(1) 



In order to take into account the limited resolving capac- 
ity of the instrumental observation facilities and to pass 
to the interval description of the system dynamics (1), 
we introduce following to [1] quantized spaces of its vari- 
ables in the form of lattices with the following periods: 
p - for radius r, a - for angle (p, v - for impulse p r and 
d - for time t. The designated periods will characterize 
the resolving capacity of corresponding observation pro- 
cedures at the theoretical level. Thus the motion of a 
particle is described by not real numbers, but integer in- 
terval of coordinates R n = R(T n ),& n — (f>(T„). impulse 
P/j = P r (T n ) and time T n in the following form 



R n =r n + [-N,N], 
$n = <Pn + [-N ) N] ) 
P r n =p n + [-N,Nl 
T n = [nN,nN + 2N], 



(2) 



where n = 0, 1, 2, ... - the number of the temporary vari- 
able lattice site with period r = AN, and N 3> 1 - the 
interval (integer) number. 

According to (2) the state of the considered system is 
localized not in a point but in a multidimensional inter- 
val with the following absolute values of the half-width: 
ay = pN - on variable r, a v — aN - on variable tp, 
f3 r = vN - on variable p r and t - on variable t. In terms 
of classical representations these values of half- width can 
be interpreted as the characterization of the numerical 
description of system. At that the given "precision" does 
not vary while the specified multidimensional interval 
moves with time in the quantized phase space forming 
an interval tube. Later we shall show that every such 
tube contains at least one classical trajectory. 

As follows from the theory [l[ , the designated values of 
the half-width are connected with the maximum veloci- 
ties of change of the system generalized coordinates and 
impulses by ratios 



tVt 



rW r 



(3) 



where V r > \dr/dt\ maXl V v > \dip/dt\ max , W r > 
\dp r I ^i|max- At that the equations of dynamics for the 



Pn+l = Pn 



p 2 rl 



<Pn+l = Pn + 



cr \p 2 ri - I J 



I = I + S^ sin(er</? m ), 



m=0 



w = W r r/V r . 

Here symbol [x] [J, 0] means the procedure of a round-off 
of real number x up to the nearest integer, and function / 
where 1$ and S are the constants defined by the problem 
intervalization parameters, means the inertia moment of 
an interval particle. The occurrence of I in the equa- 
tion for ip n of the system (4) is caused by the fact that 
the interval particles have spin, that, in its turn, and as 
is mentioned above, is a consequence of finiteness of the 
localization area size of the interval particles. The de- 
tailed discussion of this effect is beyond the problem of 
this work. It will be an objective of a special paper. 

Let's apply the equation system (4) to calculate the 
motion of the particle over the closed orbit taking for 
definiteness fi = l,c = l,e = 0.3 (e - eccentricity). 
Perform the calculations with step t = 0.07238 that 
is ~ 0.01 of the motion period over the selected orbit. 
Take N = 10 7 ,V r = 0.303, W r = 0.52, V v = 1.71, 1 = 
1.3 x 10~ 4 ,5 = 1.5 x 10~ 5 for the specified parameters. 
In this connection according to (3) the half-width of the 
interval of particle localization in the space shall not ex- 
ceed a r « 0.022. 

The results of calculations are presented in Figures 1-3. 
Fig. 1 shows the performance of full energy conserva- 
tion laws and the kinetic moment of system. The figure 
presents the relative fluctuations of current energy val- 
ues AE/\E\ and kinetic moment Ac/c for time interval 
of 2.4 x 10 6 periods. These fluctuations are computed 
with the following formulas: 

AE = E int - E, 



where 



^int - o 



4nt 



2 2 i 
V Pn + 



p 2 rl 



pr n ' 
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TABLE I: The time dependence of 
phase variables r n and p n on a pe- 
riodic trajectory. 

n r n pn_ 

350790178 -20 

1 352000368 9742590 

2 354394531 19274116 



0,004 



100 350790178 -20 

101 352000368 9742590 

102 354394531 19274116 




1x10 4 2x10" 2,4x10 s 



2,4x10° 



c int 



2 2 

P 

p 2 r 2 - 



E 



2c 2 



It follows from Fig. 1 that the designated fluctuations 
lay in strictly fixed intervals which do not vary with the 
course of time, and equal by the width to ~ 0.008 - for 
energy, and to ~ 10 -4 - for kinetic moment. 

The Fig. 2 shows effect of an interval tube. Fig. 2 
presents the time dependence of Ax = — ^kep and 
Ay = £/i n t — Z/kep characterizing the deviations of co- 
ordinates of the interval tube axis X[ n ^ — pr n cos(aip n ), 
ywjt = pr n sin(cr<^ rl ) on the corresponding values of Ke- 
pler's coordinates 



^kep 



2/kep 



c 2 cos(aip n ) 
/i He cos(a(p n ) 

c 2 sm(aip n ) 
/i 1 + e cos(aipn) 



According to Fig. 2 these deviations do not fall outside 
the limits of the fixed intervals with the half- width a x sa 
0.006 and a y sa 0.012, i.e. the specified values satisfy 
the condition of localization a x . y < a r . This implies 
that the classical (Kepler's) trajectory all over does not 
fall outside the limits of the intervals predicted by the 
equations (4), and entirely lies inside the corresponding 
interval tube. 

One more interesting property of the interval tube be- 
coming apparent in case of periodic motions is its clos- 
ability. In the example under analysis such closing is 
observed on a phase plane r n ,p n . As calculations show, 
for initial conditions 



r Q = 350790178, p Q = -20 



(5) 



the closing of the interval tube R n , (interval analogue 
of the classical trajectory r,p r ) occurs in 100 steps, i.e. in 
one turn of a particle. This result is illustrated in Table 
1 where the values of quantum numbers r n and p n are 
calculated by means of (4) in the beginning of the motion 
and in 100 steps. 

As is obvious, coincidence of these values represents 
not an approximate result, but an absolutely exact nu- 
merical result. Such property is a distinctive feature of 



FIG. 1: The time dependence between the relative quantities 
of full energy fluctuations (AE/\E\) and the orbital kinetic 
moment (Ac/c) of Kepler's particle. 




-0,009 
-0,012 



1x10 4 2x10" 2,4x10 s 
n 



2,4x10 



FIG. 2: The time dependence between the differences of the 
interval and classical values of Cartesian coordinates of Ke- 
pler's particle. 



the motion interval description and in principle is not 
realized within the limits of classical calculation means. 

The picture of the interval description of particle mo- 
tion in Kepler's Problem is not complete if not to con- 
cern one more aspect of the interval theory, namely, an 
estimation of particle position predictability horizon. As 
follows from [1|, occurrence of such characteristic as "life- 
time" of the theoretical prediction in the interval theory 
is a direct consequence of the explicit accounting of the 
interval nature of physical system phase space. The ele- 
mentary cell of such space contains not one but 2 2s points 
in the form of combinations of integer values of impulses 
and coordinates, the single-type variables in which can 
differ from each other by not more than one quantum. 
The beam 2 2s of interval tubes coming out of such cell 
(physical point) gives the reliable prediction of the cur- 
rent system state till the integration of intervals of all 
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phase variables for each of variables keeps in the interval 
of width AN. At violation of this condition if only for one 
of the variables it is possible to speak about an output 
beyond the horizon of the system state predictability. 

This is the formality of the problem. In order to char- 
acterize the problem practically, let us consider particu- 
lar procedures carried out when estimating the horizon 
of Kepler's particles position predictability. And by the 
example of these procedures we shall demonstrate tech- 
nology. 

Let's previously notice that the phase space of the 
considered system represents section p v — c and forms 
a three-dimensional lattice with periods p, v and a for 
phase variables r n ,p n and (p n , accordingly. Hence, the el- 
ementary cell of such space is a cube with a one-quantum 
side. Let's introduce an estimation of divergence of a 
beam of the trajectories coming out of this cell during 
the initial moment of time. We analyze all the possi- 
ble combinations of trajectory pairs, and select that very 
pair in which the trajectories diverge with the maximum 
velocity. In this case we calculate the divergence time 
comparing coordinates x m {. ^, y m ^ j and 2: J/Jnt 2 °^ 
each pair applying criterion 



Al < I 



loc 



(G) 



where 



Al = \j Oint,l - ^int^) 2 + faint,! ~ ^int^) 5 



and Zj oc ?s 0.024. This value of width of the interval of 
the spatial localization of the particle is selected on the 
basis of the computation results presented in Fig. 2. 

As the first example consider the dynamics of Al for 
a beam of trajectory tubes coming out of the cell, one 
top of which corresponds to the initial condition (5). For 
this condition the result of computation of the time de- 
pendence of value Al is shown in Fig. 3a and corresponds 
to the pair 



r i = 350790178, p m 
r 02 = 35 7 90 1 79, p Q2 



-20, <p 01 



-20, 



31- (?) 

-19, <p 02 = -20. 



It follows from Fig. 3a, because of (6), that the particle 
position predictability horizon is in the order of 2.5 x 10 5 
periods of the orbit motion. 

Estimating this result, it is well to bear in mind that 
it corresponds to the localization of the initial particle 
state in terms of a phase space elementary cell. In prac- 
tice such high degree of localization is not always achiev- 
able. Taking into consideration the real possibilities of 
observation tools, we should expand the area of the sys- 
tem initial state localization to some set of elementary 
cells .[T3| Such roughening of the system description leads 
as a rule to reduction of the dynamics reliable prediction 
time. This effect is shown in Figs. 3b and 3c. At that 
Fig. 3b describes pair 



r i = 350790168, p i = -20, <p 0l = -20, 
r 02 = 35 7 90 1 69, p Q2 = -19, <p 02 = -20, 



(8) 
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FIG. 3: The time dependence of Al distance between Ke- 
pler's particle locations belonging to two phase trajectories 
physically imperceptible in the initial moment. 



and Fig. 3c describes pair 

r i = 350790228, p 01 = -30, <p 01 = -20, 
ro2 = 350790227, p 02 = -29, tp 02 = -20. 



(9) 



The elementary cell of pair (8) is ten (10) quanta apart 
from the cell of pair (7) over r„ . This cell is characterized 
by an essentially smaller predictability horizon (if com- 
pare with (7)) equal to ~ 7 x 10 3 periods of the particle. 
We have still smaller predictability horizon of ~ 4.5 x 10 3 
periods for pair (9) being 50 quanta apart from pair (7) 
over r n and 10 quanta apart - over p n . 



III. ORBITAL DYNAMICS OF PLANETS OF 
SOLAR SYSTEM 

With the help of the interval equations [H numerical 
computation of motion of the major planets of the Solar 
system in the time interval of 6 x 10 9 years has been made. 
Computations have been made with step r = 0.012 year 
for N = 10 9 . The masses of the planets have been taken 
from the system of constants IAU 1964. The initial rect- 
angular heliocentric coordinates and velocities are related 
to the equator and correspond to the stage of 1949, Dec. 
30.0 ET=JED 2433280.5. To study a long-term behav- 
ior of the planet orbits and their possible drift in a chaos 
zone, the time dependence of maximum eccentricities and 
inclinations [ll| of orbits for the interval centers of the 
mentioned variables has been calculated. The values of 
the half-widths of these intervals (a e - for eccentricity 
and ai - for inclination) are presented in Table 2. The 
maximum values have been selected from the set of val- 
ues in time intervals of 6 million years, i.e. the technique 
similar to 0] has been applied. The results of the calcu- 
lations are presented in Figures 4-7. 



5 



TABLE II: The values of half- 
widths of intervals of planetary or- 
bit eccentricity (a e ) and inclination 

(«0- 

a e on 
Mercury 0.17 
Venus 6.5E-3 



Earth 

Mars 

Jupiter 

Saturn 

Uranus 



0.0165 

0.025 

0.035 

3.5E-3 

3.0E-3 



<5.0E-5 



Neptune 4.5E-3 
Pluto 0.045 



0,70 
0,65 
0,60 
0,55 

S 0,48 

| 0,46 

2 0,44 

E 0,42 

'-— 1 



£ 0,46 

S 0,45 

<g 0,44 

c 0,43 



0,45 



2 3 4 

Time (billion years) 




0,01680 
0,01674 
0,01668 



0,0936 
0,0930 
0,0924 



2 3 4 

Time (billion years) 



FIG. 4: The time dependence of maximum orbit eccentricities 
of internal planets (for interval centers). 



FIG. 5: The time dependence of maximum orbit inclinations 
of internal planets (for interval centers). 



0,059 
0,058 



-y/ 



0,0 0,1 0,2 0,3 0,4 0,5 0,6 5,4 5,5 5,6 5,7 5,8 5,9 6,0 



0,088 
0,087 
0,086 



w^^ 1 ^ (b) 



0,0 0,1 0,2 0,3 0,4 0,5 0,6 5,4 5,5 5,6 5,7 5,8 5,9 6,0 



o 0,180 
E 0,177 



0,035 
0,030 



(c) 



^Arr^^ (d) 



0,0 0,1 0,2 0,3 0,4 0,5 0,6 5,4 5,5 5,6 5,7 5,8 5,9 6,0 



(e) 



0,0 0,1 0,2 0,3 0,4 0,5 0,6 5,4 5,5 5,6 5,7 5,8 5,9 6,0 
Time (billion years) 



Figs. 4 and 5 contain graphs of the above-mentioned 
dependences for the internal planets. They show a pos- 
sible drift of orbit element values. Therewith the growth 
of the Mercury eccentricity (for which the zone of chaos 
is maximum) is limited (according to Fig. 4 and Table 
2) by an interval of values the upper limit of which does 
not exceed ^0.38. For the considered time interval this 
result agrees with the results obtained in 0, d, 0| • 

Computation of behavior of the external planets (Figs. 
6 and 7) also coincides with that having been obtained 
earlier by the authors having been mentioned and also 
by |6J|. The Uranus is an exception. For the Uranus, as 
seen from Fig. 6c, value e mas (with regard to the interval 
correction of Table 2) comes up to value of ~0.178. Thus 
Smin ~ 2 x 10~ 2 and as computations show, the Uranus 
eccentricity drift in this range has a random (stochastic) 
character, i.e. it can be regarded as a feature of chaotiza- 
tion of its motions and a motion of the external planets as 
a whole. A source of such chaotization is overlapping of 
the components of triple resonance of average motions of 



FIG. 6: The time dependence of maximum orbit eccentricities 
of external planets (for interval centers) : (a) - Jupiter, (b) - 
Saturn, (c) - Uranus, (d) - Neptune, (e) - Pluto. 



the Jupiter, Saturn and Uranus having been analyzed in 
the work of [l(| ■ Other source is overlapping of resonant 
areas in the vicinity of the Uranus and Neptune orbits 
analyzed in flD . [l~2T | . At the same time, the obtained re- 
sult is still not sufficient for final conclusions. To get an 
unambiguous answer about the nature of motion of the 
external planets, additional computational investigations 
are necessary. 

Completing the description of the results of this part of 
work, let us touch upon the problem of labor coefficient 
of the conducted computations. They have been carried 
out with the help of a personal computer with processor 
AMD Athlon 64 3400+ . To improve the reliability, the 
computations have been conducted with extended preci- 
sion having demanded the doubling of computation time. 
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FIG. 7: The time dependence of maximum orbit inclinations 
of external planets (for interval centers) : (a) - Jupiter, (b) - 
Saturn, (c) - Uranus, (d) - Neptune, (e) - Pluto. 



This time has been 1000 hours for time interval of 6 bil- 
lion years (5 x 10 12 steps). 

IV. CONCLUSIONS 

The conducted computations show the efficiency of the 
interval approach, its adequacy to the problems of ce- 



lestial bodies' dynamics. The interval means of motion 
description assure to obtain a solution in the defined 
strictly fixed interval of divergence from the classical tra- 
jectory (in case of periodic motions) and to attribute a 
unique property of absolutely exact closing. In solving 
the problems of long-term motion prediction the interval 
approach has one more advantage. Being free of round- 
off error accumulation effect it enables to computation- 
ally investigate the dynamics of planets in an arbitrary 
large time interval. At that the interval theory includes 
special computational procedure for estimating the hori- 
zon of the motion investigated aspects predictability. 



The examples considered in the paper do not obviously 
solve all the problems related to the interval planetary 
dynamics. Along with academic problems, a set of prob- 
lems being solved thanks to the creation of highly-precise 
ephemerides of planets and satellites is particularly ac- 
tual. Moreover, the potential of the interval approach 
can be realized with the greatest efficiency in the class of 
applied problems. 



It is obvious, that such realization is impossible with- 
out interest on the part of experts belonging to the ap- 
propriate application areas. And one of the problems of 
the present publication consists in turning the experts' 
attention to perspectiveness of application of the interval 
theory methods and algorithms. 
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